Electronic, Mechanical, and Piezoelectric Properties of ZnO Nanowires 
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Hexagonal [0001] nonpassivated ZnO nanowires are studied with density functional calculations. 
The band gap and Young's modulus in nanowires which are larger than those in bulk ZnO increase 
along with the decrease of the radius of nanowires. We find ZnO nanowires have larger effective 
piezoelectric constant than bulk ZnO due to their free boundary. In addition, the effective piezo- 
electric constant in small ZnO nanowires doesn't depend monotonously on the radius due to two 
competitive effects: elongation of the nanowires and increase of the ratio of surface atoms. 

PACS numbers: 77.65.-j,62.25.+g,73.22.-f,61.46.-w 



ZnOQ is one of the most important materials 
due to its three key advantages: semiconducting 
with a direct wide band gap of 3.37 eV and a 
large excitation binding energy (60 meV), piezoelec- 
tric due to non-central symmetry in the wurtzite 
structure, and biocompatible. Recently, a diversity 
group of ZnO nanostructures including nanowires 3, 
nanobeltsi^, nanospringS|^, nanocombs 5j, nanorings 3j 
nanobowsQ, and nanohelices^, 0| have been synthe- 
sized under specific growth conditions. ZnO nanos- 
tructures could have novel applications due to their 
unique physical and chemical properties arising from 
surface and quantum confinement. In particular, ZnO 
nanowires with relatively simple structures are impor- 
tant one-dimensional (ID) nanostructures. Experimen- 
tally, the group of Wang had synthesized well-aligned 
[0001] ZnO nanowires enclosed by facet {1010} surfaces 
[lol fllj . Room-temperature ultraviolet lasingjl2| and 
piezoelectric nanogenerators based on ZnO nanowire ar- 
rays have been demonstrated^^]- Rectifying diodes of 
single ZnO nanobelt/nanowire-based devices and a 
ZnO nanowire photodetectorfisj were fabricated very re- 
cently. 

Although many studies on ZnO nanowires have been 
conducted, there are some important issues remained to 
be addressed. First, the mechanical properties, especially 
the Young's modulus of ZnO nanowires are on debate in 
the literature© lii llM llSl E| • 

For instance, Chen et al. 
showed that the Young' modulus of ZnO nanowire 
with diameters smaller than about 120 nm is signifi- 
cantly higher than that of bulk ZnO. However, the elas- 
tic modulus of vertically aligned [0001] ZnO nanowires 
with an average diameter of 45 nm measured by atomic 
force microsc opy was found to be far smaller than that 
of bulk ZnO 1171. The second issue is about the elec- 
tromechanical coupling in ZnO nanowires. The effective 
piezoelectric coefficient of individual (0001) surface dom- 
inated ZnO nanobelts measured by piezoresponse force 
microscopy was found to be much larger than the value 
for bulk wurtzite ZnO ^3]. In contrast. Fan et al. showed 
that the piezoelectric coefficient for ZnO nanopillar with 



the diameter about 300 nm is smaller than the bulk 
values 22]. They suggested that the reduced electrome- 
chanical response might be due to structural defects in 
the pillars [2^ • Whether the electromechanical coupling 
is enhanced or depressed in defect-free ZnO nanowires 
is not clear. Thirdly, although it is well known that the 
quantum confinement effect will decrease the band gap of 
passivated nanowires, the question that how the dangling 
bond in bare ZnO nanowires affects the band gap remains 
open. The fundamental study on these issues is crucial 
for developing future applications of ZnO nanowires. 

In this letter, we have studied the electronic, mechani- 
cal, and piezoelectric properties of [0001] ZnO nanowires 
using first-principles methods for the first time. We find 
that the band gap increases along with the decrease of 
the radius of ZnO nanowires due to the radial confine- 
ment. The Young's modulus of nanowires is larger than 
bulk ZnO, in agreement with the experimental results of 
Chen et alJld\. The effective piezoelectric constant in 
ZnO nanowires is larger than that of bulk ZnO due to 
the free boundary of nanowires. Moreover, the effective 
piezoelectric constant in small ZnO nanowires doesn't de- 
pend monotonously on the radius due to two competitive 
effects. 

Our calculations are performed using the SIESTA 
package [2^. a standard Kohn-Sham density- functional 
program using norm- conserving pseudopotentials and nu- 
merical atomic orbitals as basis sets. The local den- 
sity approximation (LDA) |24| to the exchange correlation 
functional is employed. The DZP basis sets for both Zn 
and O with the energy shift parameter 20 mRy are used. 
The mesh cutoff parameter for real space integration is 
set to 350 Ry to obtain accurate atomic forces. Bulk 
ZnO has a wurtzite structure with the noncentral sym- 
metry, resulting in a normal dipole moment and spon- 
taneous polarization along the c-axis. The computed 
(experimental j25|) lattice parameters of bulk wurtzite 
ZnO are a = 3.17 A(3.25 A), c = 5.18 (5.20 A), u = 0.374 
(0.381). The calculated LDA direct band gap is 0.63 eV, 
in good agreement with other LDA results j26l 1271]. 

Here we mainly focus on [0001] ZnO nanowires 
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TABLE I: Unrelaxed diameter (Do) relaxed diameter (D), 
relaxed lattice constant (c), strain energy (_E), band gap (-Eg), 
Young's modulus (-E3), and effective piezoelectric constant 
(633) of five different ZnO nanowires ( A, B, C, D, and E) 
and bulk ZnO. 



A 



B 



C 



D 



E bulk 



Do (A) 3.66 9.68 15.96 22.27 28.59 

D (A) 3.32 9.32 15.61 21.97 28.33 

c (A) 5.335 5.302 5.270 5.234 5.215 5.180 

Eg (eV) 2.40 1.54 1.09 0.85 0.75 0.63 

E3 (GPa) 363 242 217 189 182 147 

6^3 (10~^VcA/ion) 2025 1837 1879 1986 1961 1453 



since almost all ZnO nanowires grown along the [0001] 
directionfiollllt. The nanowires are enclosed by six facet 
{1010} surfaces with low surface energy. In wurtzite 
bulk ZnO, each Zn (O) bonds with four O (Zn) atoms. 
In nanowires, the surface Zn (O) atom is only bound 
to three nearest neighbours. In this study, we consider 
only bare ZnO nanowires without saturating the surface 
atoms. The nanowires are modeled by hexagonal super- 
cells whose lateral lattice constants are so large that there 
is almost no interaction between the nanowires. Five 
ZnO nanowires with diameter (Here the diameter is de- 
fined as the largest lateral distance between atoms) rang- 
ing from about 0.3 to 2.8 nm, labeled as A, B, C, D, and 
E respectively, are examined. The largest nanowire E 
contains 300 atoms in the unitcell. Since the lattice con- 
stant of small ZnO nanowires might differ significantly 
from the bulk counterpart, we optimize both the lattice 
constant c and internal coordinates of nanowires. The 
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FIG. 1: (Color online) Relaxed structures for ZnO nanowires 
with different radius. We label these nanowires as A, B, C, 
D, and E. 



relaxed structures of these ZnO nanowires are shown in 
Fig. ^ Since both surface Zn and O atoms move inwards 
and Zn atoms move much more, it looks like surface O 
atoms rotate outwards. The diameters of relaxed and 
unrelaxed nanowires are shown in Tabled We can see 
that the diameters of relaxed nanowires are smaller than 
those of unrelaxed nanowires by almost 0.3 A. The re- 
laxation of surface atoms in ZnO nanowires is similar to 
that in ZnO [1010] surface [isj]. Along with the shrinkage 
of surface atoms, the lattice contant c of ZnO nanowires 
is increased when compared with that of bulk ZnO, as 
are shown in Tabled We can see that the elongation of 
small nanowires is considerably large, however, the lat- 
tice constants of large nanowires tend to approach that 
of bulk ZnO. 

We have calculated the electronic structures of these 
ZnO nanowires. All [0001] ZnO nanowires are found to 
be semiconducting. Although the LDA usually under- 
estimates the band gap, the trend of the band gaps of 
ZnO nanowires predicted from the LDA calculations are 
expected to be correct. The LDA band gaps are shown 
in Table |l| Clearly, the band gaps of ZnO nanowires 
increase monotonously along with the decrease of the ra- 
dius of nanowires. In comparison with bulk ZnO, the 
increment of the band gap of nanowire A can be as large 
as 1.77 eV. The blueshift of the band gap should be due 
to the quantum confinement effect. However, the gap- 
broadening effect in nonpassivated nanowires is unusual. 
For example, nonpassivated Si nanowires grown along the 
(100) direction are found to be metallic and semimetallic 
due to the presence of surface states [i^. To see the de- 
tails of the electronic structure of ZnO nanowires, we plot 
the band structure of nanowire B in Fig.[2Ia). The band 
structures of other ZnO nanowires are similar to that of 
nanowire B. Clearly, ZnO nanowires have a direct gap 
at F. To gain an insight into the character of the lowest 
unoccupied molecular orbital (LUMO) and highest oc- 
cupied molecular orbital (HOMO), we show the charge 
density of the LUMO and HOMO states at F in Fig.l^Ib) 
and (c) , respectively. We can see that the LUMO is delo- 
calized in the whole nanowire, indicating that it is a bulk 
state instead of a surface state. The delocalized distri- 
bution is also responsible for the large dispersion of the 
LUMO from F to A. Furthermore, due to the delocalized 
character, the LUMO energy will increase substantially 
in small radius nanowires due to the radial confinement. 
Since the charge density of the HOMO at F mainly con- 
tributed by surface atoms, the HOMO is a surface state. 
The HOMO at F lies only 80 meV above the top of va- 
lence band of bulk ZnO, and its position changes little 
in nanowires with different diameters since the HOMO 
is mainly composed by surface O 2p like dangling bonds. 
The HOMO is similar to the PI surface state in the ZnO 
(1010) surface as identified by Schroer et al. 28]. It is 
the different response to the radial confinement between 
the LUMO and HOMO that leads to the overall gap- 
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broadening effect in ZnO nanowires. 




FIG. 2: (Color online) (a) The band structure of nanowire 
B. (b) The charge density of the LUMO state at F for 
nanowire B. (c) The charge density of the HOMO state at 
r for nanowire B. 

Now wc turn to study the mechanical properties of 
ZnO nanowires. An important mechanical parameter for 
describing a one-dimensional system is the Young's mod- 
ulus. We calculate in the [0001] direction by using 
the following formula: 

E, ^ ^<fElde\ (1) 

where E is the total energy, e is the axial strain, and V 
is the system volume. For nanowires, we define the vol- 
ume as V = S X c, where S is the area of the cross 
section of the nanowire and c is the lattice constant. 
For each nanowire, we perform seven calculations with 
e range from —3% to 3%. In each calculations for the 
nanowires, we fully relaxed the internal coordinates. To 
compare our results for nanowires with bulk ZnO, we 
also calculate the Young's modulus of bulk ZnO. In the 
calculations, only the lattice constant c is fixed to the 
strained value, all other degrees of freedom including 
the lateral lattice constants are optimized. Our results 
for bulk ZnO and ZnO nanowires are reported in Ta- 
ble n The calculated Young's modulus of bulk ZnO is 
147 GPa, which is close to the value 140 GPa deduced 
from the experimental elastic constants!^. For ZnO 
nanowires, the Young's modulus increase monotonously 
along with the decrease of the radius of nanowires and are 
larger than that of bulk ZnO. To estimate the Young's 
modulus of larger ZnO nanowires, the calculated re- 
sults are fitted using Es{D) = 145.87 1172.63D-1 - 
1477.58I?-2 -I- 4830.55LI-3, where Young's modulus E-i 
is in unit of GPa and diameter D is in unit of A. The 
trend of Young's modulus accords qualitatively with the 



experimental results[l^. Our results also agree qualita- 
tively with the Young's modulus of ZnO nanobelts with 
rectangular cross-sections in the [0001] orientation calcu- 
lated from empirical molecular dynamics simulations |2nj]. 
However, the experimental Young's modulus of ZnO 
nanowires are significantly larger than our results. For 
example, the Young's modulus of a ZnO nanowire with 
the diameter 20 nm estimated using their fitted formula 
is 202 Gpa, however, our value estimated from the fitted 
formula is only 152 Gpa. The reason for the discrepancy 
is not very clear. 

Before we proceed to discuss the piezoelectric effect 
in ZnO nanowires, we first calculate piezoelectric con- 
stant 633 of bulk ZnO. The piezoelectric constants e^j is 
defined as follows: — dVijdtj^ where P is the to- 
tal polarization, and tj is the strain tensor component. 
Here, the piezoelectric constants are computed by com- 
bining the Berry phase method in the modern theory of 
polarization|3lj with the finite difference method. The 
piezoelectric effect in tetrahedrally bonded semiconduc- 
tors results from two different terms of opposite sign: The 
"clamped-ion" and the "internal-strain" contributions. 
The calculated clamped-ion and relaxed-ion piezoelectric 
constant 633 for bulk ZnO is —0.77 and 1.29 C/m^ re- 
spectively, which accords with others' results (—0.75 and 
1.28 C/m^)!!!. 

As the conventional definition of piezoelectric constant 
for three-dimensional bulk is not appropriate for describ- 
ing the piezoelectric properties of one-dimensional sys- 
tems, we define the atomic averaged effective piezoelec- 
tric constant as: 633=633 x ^gcell/-^' where N is the 
number of atoms, and V^cgH is the volume of the super- 
cell. The results are presented in Fig. Eta) and Table 
We can clearly see that 633 of nanowires is considerably 
larger than that of bulk ZnO. We find that larger 633 in 
nanowires is caused by the free boundary of nanowires. 
In the calculations of 633 , the lateral lattice constants are 
fixed. When bulk ZnO is strained along the c axis, atoms 
can not relaxed freely along the lateral directions. How- 
ever, when nanowires are compressed or elongated along 
the c axis, atoms can relaxed freely along the lateral di- 
rections due to the free boundary. Hence, the effective 
piezoelectric constant 6^3 in nanowires should approach 
633 — 26^]^ X where 633 and A^-^ is effective piezoelectric 
constants of bulk ZnO, v — —61/63 — —62/63 is Poisson's 
ratio (61 = 62 = (a — ao)/ao, 63 = (c — co)/co, a is the 
relaxed lateral lattice constant when the lattice constant 
of ZnO is changed to c). Since > and 63^ < 2(|, 
633 in nanowires is larger than 6^3. For the same rea- 
son, effective piezoelectric coefficient in ZnO nanobelts 
was also found experimentally to be larger than that of 
bulk ZnO by Zhao el aL^|. One might expect smaller 
nanowires have larger piezoelectric constant due to their 
larger ratio of surface atoms. However, Fig. El^a) shows 
that 633 doesn't monotonously depend on the radius of 
ZnO nanowires: 633 for nanowire D is larger than that 
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for nanowire C and E. We attribute this anormal behav- 
ior to the increase of the lattice constant along with the 
decrease of the radius of ZnO nanowires, since the piezo- 
electric constant of a certain nanowire decreases along 
with the increase of the lattice constant c, as shown in 
Fig-EJb). In fact, we also calculate the effective piezo- 
electric constant of the nanowires with bulk lattice con- 
stant c, and the results (also shown in Fig.O^a)) indicate 
a monotonous decreasing dependence of 633 upon the ra- 
dius of nanowires. From Fig. |2Ia), we can also see that 
the difference of 633 between nanowires D and E with the 
lattice constant c relaxed or fixed is very small, suggest- 
ing that larger nanowires will have similar 633 as nanowire 
E. 
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FIG. 3: (a) The dependence of the effective piezoelectric con- 
stant 633 upon the diameters of the nanowires with optimized 
lattice constant c and nanowires with bulk lattice constant c. 
The horizontal dotted line indicates the piezoelectric constant 
of bulk ZnO. (b) shows the polarization of ZnO nanowire B 
with different c. 

To summarize, we have carried out comprehensive 
first-principles studies on [0001] ZnO nanowires. The 
band gaps of ZnO nanowires increase along with the 
decrease of the radius of nanowires due to the quan- 
tum confinement effect. The Young's modulus of thiner 
nanowires is larger than that of thicker nanowires. Our 
calculations indicate that the effective piezoelectric con- 
stant 6^3 of ZnO nanowires is larger than that of bulk 
ZnO. In addition, wc find a non-trivial dependence of 
the electromechanical coupling of ZnO nanowires upon 
the radius as a result of the competition between two 
opposite factors. Our results support the application of 
ZnO nanowires as nanosensors and nanoactuators. 

This work is partially supported by the National Nat- 
ural Science Foundation of China (50121202, 20533030, 
10474087), by the USTC-HP HFC project, and by the 
SCCAS and Shanghai Supercomputer Center. 



* Corresponding author. E-mail: jlyang@ustc.edu.cn 
[1] U. Ozgiir, Ya. I. Alivov, C. Liu, A. Teke, M. A. 
Reshchikov, S. Dogan, V. Avrutin, S. J. Cho, and H. 



[2: 

[s: 

[4: 
[s: 

[6: 

[t: 

[s: 

[9: 

[10: 

[11 

[12 

[13 
[14 

[15 

[16 

[17 

[18 

[19 

[20 

[21 

[22; 

[23; 

[24; 

[25; 

[26; 

[27 

[28; 

[29 

[30 
[31 



Morkoc, J. Appl. Phys. 98, 041301 (2005). 

M. H. Huang, Y. Wu, H. Feick, N. Tran, E. Weber, and 

P. D. Yang, Adv. Mater. 13, 113 (2001). 

Z. W. Pan, Z. R. Dai, and Z. L. Wang, Science 291, 1947 

(2001). 

X. Y. Kong and Z. L. Wang, Nano Lett. 3, 1625 (2003). 
Z. L. Wang, X. Y. Kong, and J. M. Zuo, Phys. Rev. Lett. 
91, 185502 (2003). 

X. Y. Kong, Y. Ding, R. S. Yang, and Z. L. Wang, Science 
303, 1348 (2004). 

W. Hughes and Z. L. Wang, J. Am. Chem. Soc. 126, 
2709 (2004). 

R. S. Yang, Y. Ding, and Z. L. Wang, Nano Lett. 4,1309 

(2004) . 

P. M. Gao, Y. Ding, W. J. Mai, W. L. Hughes, C. S. Lao, 

and Z. L. Wang, Science 309, 1700 (2005). 

P. X. Gao, Y. Ding, and Z. L. Wang, Nano Lett. 3, 1315 

(2003). 

X. D. Wang, C. J. Summers, and Z. L. Wang, Nano Lett. 

4, 423 (2004). 

M. H. Huang, S. Mao, H. Feick, H. Q. Yan, Y. Y. Wu, 
H. Kind, E. Weber, R. Russo, and P. D. Yang, Science 
292, 1897 (2001). 

Z. L. Wang and J. Song, Science 312, 242 (2006). 
C. S. Lao, J. Liu, P. X. Gao, L. Y. Zhang, D. Davidovic, 
R. Tummala, and Z. L. Wang, Nano Lett. 6, 263 (2006). 
J. B. K. Law and J. T. L. Thong, Appl. Phys. Lett. 88, 
133114 (2006). 

C. Q. Chen, Y. Shi, Y. S. Zhang, J. Zhu, and Y. J. Yan, 
Phys. Rev. Lett. 96, 075505 (2006). 

X. D. Bai, P. X. Gao, Z. L. Wang, and E. G. Wang, Appl. 
Phys. Lett. 82, 4806 (2003). 

K. Yum, Z. Wang, A. P. Suryavanshi, and M. F. Yu, J. 
Appl. Phys. 96, 3933 (2004). 

J. Song, X. Wang, E. Riedo, and Z. L. Wang, Nano Lett. 

5, 1954 (2005). 

A. J. Kulkarni, M. Zhou, and F. J. Ke, Nanotechnology 
16, 2749 (2005). 

M. H. Zhao, Z. L. Wang, and S. X. Mao, Nano Lett. 4, 
587(2004). 

H. J. Fan, W. Lee, R. Hauschild, M. Alexe, G. L. Rhun, 
R. Scholz, A. Dadgar, K. Nielsch, H. Kalt, A. Krost, M. 
Zacharias, and U. Gosele, SmaU 2, 561 (2006). 

J. M. Soler, E. Artacho, J. D. Gale, A. Garcia, J. Jun- 
quera, P. Ordejon, and D. Sanchez-Portal, J. Phys.: Con- 
dens. Matter 14, 2745 (2002). 

D. M. Geperley and B. J. Alder, Phys. Rev. Lett. 45, 
566(1980); J. P. Perdew and A. Zunger, Phys. Rev. B 
23, 5048 (1981). 

K. Yoshio, A. Onodera, H. Satoh, N. Sakagami, and H. 
Yamashita, Ferroelectrics 264, 133 (2001). 
X. Wu, D. Vanderbilt, and D. R. Hamann, Phys. Rev. B 
72, 035105 (2005). 

M. Usuda, N. Hamada, T. Kotani, and M. van Schilf- 
gaarde, Phys. Rev. B 66, 125101 (2002). 
P. Schroer, P. Kriiger, and J. PoUmann, Phys. Rev. B 
49, 17092 (1994). 

R. Rurah and N. Lorente, Phys. Rev. Lett. 94, 026805 

(2005) . 

I. B. Kobiakov, Solid State Commun. 35, 305 (1980). 
R. D. King-Smith and D. Vanderbih, Phys. Rev. B 47, 
1651 (1993); R. Resta, Rev. Mod. Phys. 66, 899 (1994). 



